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Abstract 


The production cross sections of the B+, and Bg mesons, and of their charge con¬ 
jugates, are measured via exclusive hadronic decays in pPb collisions at the center-of- 
mass energy = 5.02 TeV with the CMS detector at the CERN LHC. The dataset 
used for this analysis corresponds to an integrated luminosity of 34.6nb^^. The pro¬ 
duction cross sections are measured in the transverse momentum range between 10 
and 60GeV/c. No significant modification is observed compared to proton-proton 
perturbative QCD calculations scaled by the number of incoherent nucleon-nucleon 
collisions. These results provide a baseline for the study of in-medium b quark energy 
loss in PbPb collisions. 
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Relativistic heavy ion collisions allow the study of quantum chromod}mamics (QCD) at very 
high temperature and density. Under such extreme conditions, a strongly interacting state 
consisting of deconfined quarks and gluons, the quark-gluon plasma (QGP) lUJEl, is predicted 
by lattice QCD calculations El- Hard-scattered partons are expected to lose energy as they 
traverse the QGP via elastic collisions and medium-induced gluon radiation. The resulting 
reduction of the measured yield of hadrons, compared to expectations based on proton-proton 
(pp) data, is often referred to as "jet quenching" ||ll|5l. The flavor dependence of jet quenching 
is one of the most important testing grounds for energy loss models Il6l- [l0l . However, other 
phenomena can affect the yield of heavy-flavor particles, independently of the presence of a 
deconfined partonic medium. For instance, modifications of the parton distribution functions 
(PDFs) in the nucleus with respect to nucleon PDFs IfTlT - IT^ could change the production rate. 
Therefore, a complete understanding of the interactions of heavy quarks in the deconfined 
medium formed in heavy ion collisions requires a thorough knowledge of their production in 
proton- (or deuteron-) nucleus, p(d)A, collisions. 

Gurrently, published data for heavy flavor production in p(d)A exist for open charm both at 
RHIG, in dAu collisions at = 200 GeV, and at the LHG in pPb collisions at = 

5.02 TeV. At RHIG, the STAR Gollaboration measured the charm spectra in the rapidity interval 
\y\ <1 from direct reconstruction of D*^ meson and from indirect electron and positron mea¬ 
surements of charm semileptonic decays I1T4I . The measured yields were found to be consistent 
within the uncertainties with the hypothesis of binary scaling (no modification with respect to 
nucleon-nucleon (NN) cross section scaled by the number of incoherent NN binary collisions). 
However, the PHENIX Gollaboration measured a significant enhancement of the production 
of heavy-flavor decay electrons in |i/| < 0.35 in high-multiplicity dAu events with respect to 
a combined, data and theory, pp reference 1(151 . Recently, PHENIX also measured a signifi¬ 
cant enhancement of heavy-flavor production via single-muon detection at backward rapidity 
(the Au-gomg direction), and a suppression at forward rapidities (the d-going direction) 1(1^ . 
This measured difference in heavy-flavor production between forward and backward rapidi¬ 
ties is significantly larger than predicted by leading-order perturbative QGD calculations with 
nuclear PDEs IITtI . In pPb collisions at the LHG, the ALIGE Gollaboration measured the pro¬ 
duction of the D meson in the —0.96 < y < 0.04 interval, and found it to be, within uncer¬ 
tainties, compatible with pp data scaled by the number of binary NN collisions, over a large 
transverse momentum (pj) range IllSl . The LHG results are well described by theoretical cal¬ 
culations that do not require a deconfined medium to be formed in the collision. This supports 
the idea that the D meson suppression at high-px observed in PbPb collisions by the ALIGE 
Gollaboration ||19l is due to parton interactions with the deconfined medium. While measure¬ 
ments, both at RHIG and LHG, support that most of the suppression observed in AA collisions 
is due to partonic energy loss, the details of the phenomena affecting open charm in pA and 
AA collisions are still to be understood. 

The production of B mesons was studied at the LHG in proton-proton (pp) collisions at = 
7 TeV over wide px and rapidity intervals by GMS (SoHZll, ATLAS Il23ll . and LHGb UMI . In PbPb 
collisions, GMS measured the non-prompt J/ip from B hadron decays at = 2.76 TeV Il25l . 
and observed a strong suppression with respect to the hypothesis of binary scaling. In this 
Letter, we extend the study of heavy-quark production in p(d)A collisions by performing the 
first measurement of exclusive B meson decays in pPb collisions. 

The B mesons are measured in a region |i/iab| < 2.4 via the full reconstruction of their decay 
channels: B+ —> ]/tp K+ —)■ p+p^K+ with branching fraction B = (6.12 ±0.19) x 10^^, B*^ — > 
]/xp K*(892) ^ y+p-K+n- with B = (5.24 ± 0.24) x 10^^ and B^ ^ ]/tp (p ^ p+p-R+R- 
with B = (3.12 ± 0.27) x 10^^ As this analysis does not separate B+ from B^, B'^ from 
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or Bg from Bg, mesons are referred to generically as B+, B^, and Bg, respectively, for the 
purposes of reconstruction. For the final cross section values, the combined results are divided 
by two to obtain an average. 

The CMS detector has excellent capabilities to reconstruct B meson decays due to the highly 
efficient muon detection system and the high-resolution silicon tracker Il27l . The data sample 
used in this analysis corresponds to an integrated luminosity of (34.6 ± 1.2) nb~^ [28l . The 
direction of the proton beam was initially opposite to the positive direction of the CMS lon¬ 
gitudinal axis Il27l , and it was reversed after 60% of the data were taken. The beam energies 
were 4TeV for protons and 1.58 TeV per nucleon for lead nuclei, resulting in a nucleon-nucleon 
center-of-mass energy of = 5.02 TeV. Because of the energy difference of the colliding 

beams, the nucleon-nucleon center-of-mass frame in pPb collisions was not at rest with respect 
to the laboratory frame. The results presented here use the convention that the proton-going 
side corresponds to positive pseudorapidity. This implies that massless particles emitted at 
pseudorapidity rjcu in the NN center-of-mass frame are detected at = rjcu + 0.465. 

A detailed description of the CMS experiment and coordinate system can be found in Ref. ETII . 
Only the detector subsystems most relevant for this analysis are described here. Charged par¬ 
ticles (tracks) are reconstructed within the range |?/iab| < 2.5 by using the silicon tracker de¬ 
tector, located in the 3.8 T magnetic field of a superconducting solenoid. Muons are identi¬ 
fied in the interval |^iab| < 2.4 with gas-ionization detectors made of three technologies: drift 
tubes, cathode strip chambers, and resistive plate chambers, embedded in the steel flux-return 
yoke of the magnet. The CMS apparatus also has extensive forward calorimetry, including two 
steel and quartz-fiber Cherenkov hadron forward (HP) calorimeters, which cover the range 
2.9 < |?/iab| < 5.2. 

Events used in the measurement are collected with a trigger requiring the presence of a muon 
with px > 3 GeV/c. To select inelastic hadronic interactions, the offline analysis requires a coin¬ 
cidence of at least one of the HP calorimeter towers (with more than 3 GeV of total energy) from 
each side of the interaction point. Events are further required to have at least one reconstructed 
primary vertex, formed by at least two tracks, with a distance from the center of the nominal 
interaction region of less than 15 cm along the beam axis. 

Several Monte Garlo (MG) simulated event samples are used to evaluate background compo¬ 
nents and signal efficiencies, specifically: i) an inclusive (prompt and non-prompt) J/ip sample; 
ii) a sample containing all B mesons decaying into a ]/xp; iii) a signal-only sample with the B+, 
B°, and Bg decays included in the present analysis. First, proton-proton collisions are simu¬ 
lated with PYTHIA 6.424 Il29ll tune Z2 Il30l and propagated through the GMS detector using the 
Geant4 package 1311 . The B meson decays are simulated with the evtgen package 1321 , and 
final state photon radiation in the B decays is simulated by PHOTOS |j33|. Then, the PYTHIA 
events are embedded into simulated pPb events produced by the HIJING generator version 
1.383 f34l , which is tuned to reproduce global event properties such as charged-hadron pj 
spectra and particle multiplicity. 

Muons are required to be within the following kinematic region: pj > 3.3GeV/cfor < 1.3, 
total momentum > 2.9GeV/c for 1.3 < < 2.2, or pj > 1.5GeV/c for 2.2 < < 

2.4 1351 . This acceptance selection is chosen so as to guarantee a single-muon detection prob¬ 
ability exceeding about 10%. Two muons of opposite charge with an invariant mass within 
150 MeV/c^ of the world-average ]/ip mass ll^ are selected to reconstruct a ]/tp candidate, with 
a mass resolution of typically 18-55 MeV/c^, degrading as a function of the dimuon rapidity. 
The B meson candidates are formed by combining ]/xp candidates with charged tracks. With¬ 
out using particle identification, assumptions need to be made about the masses of the charged 
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Figure 1: Invariant mass distributions of B+ + B (left), B° + B° (center), and B° + B° (right) 
candidates in the transverse momentum regions 10-15, 10-15, and 10-60 GeV/c, respectively. 
See the text for details. 


tracks. In calculating the mass of the B+ candidates, the single charged particle is always as¬ 
sumed to have the mass of a kaon. In the B*^ case, two invariant mass values are computed, 
corresponding to the two possible assignments of the kaon and pion masses to the two-track 
system. For B® candidates, the two charged tracks are always assumed to be kaons. Single 
track low pj thresholds of 0.9, 0.7, and 0.4GeV/c are applied in the B+, B°, and B® analyses, re¬ 
spectively, to reduce the combinatorial background, which is further minimized by additional 
selection criteria. In particular, B candidates are selected according to the probability of the 
decay vertex (the probability for the ]/tp muon tracks and the other charged track to point to 
a common vertex), the 3D flight distance (normalized by its uncertainty) between the primary 
and decay vertices, and the pointing angle, which is defined as the angle between the line con¬ 
necting the primary and decay vertices and the momentum vector of the B meson in the plane 
transverse to the beam direction. The selection is optimized for each meson species using a 
multivariate technique that uses the genetics algorithm Il36]l . in order to maximize the statisti¬ 
cal significance of the B meson signals. In the B^ and Bg analyses, the invariant masses of the 
K+tt^ and the K+K^ are required to be compatible with the masses of the K'^*(892), K*(892) 
and the (p resonances, respectively. If more than one candidate in a given event survives all the 
aforementioned selection criteria, the candidate with the best vertex x^ probability is selected. 

The raw yields of B+, B®, and Bg are extracted using a binned maximum likelihood fit to the 
B meson invariant-mass distributions in the mass range 5 < mB < 6GeV/c^. The invariant 
mass distributions of B+, B°, and Bg candidates in the pj regions 10-15,10-15, and 10-60 GeV/c, 
respectively, are shown in Fig.j^ In the case of B+ and B*^, this choice corresponds to the lowest 
Pj interval used in the analysis, while for Bg it is the only interval. The signal shape is modeled 
by two Gaussians with the same mean values (a free parameter in the fit) and different widths 
determined in MG simulations. The background is dominated by random combinations of 
prompt and non-prompt ]/xp candidates with extra particles. This combinatorial background 
is modeled by a first-order polynomial in the B+ and B'^ analyses, and by a second-order poly¬ 
nomial in the Bg analysis, as suggested by studies on the embedded inclusive ]/tp sample. The 
background component shown as a crosshatched histogram and labeled as B —> J/i/; X in Fig.j^ 
is due to mis-reconstructed B meson decays that produce broad peaking structures in the in¬ 
variant mass region below 5.4 GeV/c^. As an example, in the B+ analysis, a peaking background 
structure is created by B° —> J/t/; K*(892) decays in which one decay product is lost in the B 
candidate reconstruction. These background sources are studied with the embedded MG sam¬ 
ple including all B meson decays into final states with a ]/ip, and found to be well described 
by a superposition of 4 and 2 Gaussian functions in the B+ and B^ analyses, respectively. The 
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Figure 2: The pj-differential production cross section of B+ (left), (center), and Bg (right) 
measured in pPb collisions at = 5.02 TeV. The vertical bars (boxes) correspond to statis¬ 
tical (systematic) uncertainties. The global systematic uncertainty, listed in each panel and not 
included in the data points, comprises the uncertainties in the integrated luminosity measure¬ 
ment and the B meson branching fractions. Results are compared to FONLL calculations 1137F 
[39] , scaled by the number of binary NN collisions, represented by a continuous histogram. The 
dashed histograms represent the theoretical uncertainties for the FONLL reference. 


resulting functional form, with the overall normalization left floating, is included in the global 
fit function. This background component is found to be negligible in the Bg analysis as a con¬ 
sequence of the selection on the mass of the cp candidate. 

The px-differential production cross section of the various B meson species is computed in each 
Pt interval: 

dc^ ^ 1 1 ^^(FT)|yi3b|<2.4 

dpT |y,^,|<2.4 2ApT(Acce)|y,^^|<2.4B/:- 

^^PT)|i/[ 3 bl< 2.4 signal yield extracted in each pj interval of width Apx, ( Acce)|j^j^|^|<- 2,4 

represents the corresponding acceptance times efficiency, C is the integrated luminosity, and B 
is the branching fraction of the decay chain. The factor 1/2 accounts for the fact that the yields 
were measured for particles and antiparticles added together, but the cross section is given 
for particles only. An analogous expression holds for the rapidity-differential cross section. 
The (Acce) correction factors are evaluated using pythia+hijing simulations in each px and 
I y lab I interval, to account for the loss of signal due to the detector coverage, and to the trigger, 
reconstruction and offline selection. They vary, over the measured px and y intervals, from 
9%to 37% (3% to 15%) for B+ (B°), and they equal 8% in the single bin used for the Bg. 

The cross sections are affected by several sources of systematic uncertainties arising from the 
signal extraction, acceptance and efficiency corrections, branching fractions, and integrated 
luminosity determination. The uncertainty from the fitting procedure (varying from 10% to 
15% across all analysis intervals, for all three mesons) is evaluated by varying the probability 
distribution functions used to model the signal and background distributions. As an alterna¬ 
tive combinatorial background shape, a second-order polynomial is used for B+ and B^, and a 
third-order polynomial for Bg. The uncertainty on the signal is evaluated by considering three 
fit variations: (i) leaving free the width parameters; (ii) varying the width parameters by ±20% 
with respect to the MC value; (iii) using only one Gaussian. The maximum of all variations is 
propagated as systematic uncertainty, and it is given in the case of B+ by variation (i), and by 
variation (ii) for the other two mesons. The systematic uncertainties associated with the bin-by¬ 
bin acceptance correction (0.2% to 5.6%) are estimated by varying the shape of the generated 
B meson px and y spectra within limits defined by differences (including their statistical uncer- 
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Figure 3: The nuclear modification factors of B+ (left), (center), B® (right) mea¬ 

sured in pPb collisions at = 5.02 TeV. The statistical and systematic uncertainties on the 
pPb data are shown as bars and yellow boxes around the data points, respectively The system¬ 
atic uncertainties from the FONLL predictions are plotted separately as open blue boxes. The 
global systematic uncertainties are shown as full grey boxes at unity, and are not included in 
the data points. 


tainties) between data and MC calculations. For all three mesons, the B+ and B° px spectrum 
shapes are assumed, while only the B+ is used for the y shape. Using these shape variations, 
simplified ("toy") MC simulations are used to recalculate the acceptance in each kinematic 
bin, the maximum variation between the nominal acceptance and the toys being propagated 
as the systematic uncertainty. The systematic uncertainty due to the selection of the B meson 
candidates (4% to 11%) is equal to 1 minus the ratio of the selection efficiencies (the ratio of 
the extracted yield with and without applying the selection) estimated in data and simulation. 
In addition, an uncertainty associated with the accuracy of the best candidate selection (3%) , 
which depends on the number of reconstructed B meson candidates, is assigned. This is eval¬ 
uated by reweighting the population of the pythia+hijing events so that the distribution of 
the number of B meson candidates per event matches the one from data. The uncertainties in 
the muon trigger, and muon track reconstruction and identification efficiencies (4.5% to 7.3%), 
are evaluated by using the "tag-and-probe" technique Il40l on pPb data and the embedded MC 
sample. The systematic uncertainty associated with the track reconstruction efficiency (3.9% 
per hadronic track Bil l is estimated from a comparison of two-body and four-body de¬ 
cays in pp data and MC calculations, all samples being reconstructed with the same tracking 
algorithm as the pPb sample. The systematic uncertainty in the cross section measurement 
is computed point-by-point as the sum in quadrature of the different contributions mentioned 
above. In addition, a global systematic uncertainty is calculated to account for the uncertainties 
in the integrated luminosity value (3.5% f28|), and in the B meson branching fractions (3.1%, 
4.6%, and 8.7% for B+, B°, and Bg, respectively Il26l l. 

In Fig. |2 the px'differential production cross sections of all three B mesons measured in the 
interval] 1/iabI < 2.4 are presented, with data points placed at the center of each bin. They are 
compared to the pp cross sections obtained from fixed-order plus next-to-leading-logarithm 
(FONLL) calculations Il39l , which reproduce the B meson pj-differential cross sections in pp 
collisions at 7 TeV ||20] - IMI . The individual cross sections are obtained by scaling the FONLL 
total beauty production Il37l - l39l by the world-average production fractions of B+, B°, and Bg 
(40.2%, 40.2%, and 10.5%, respectively |[26 l). The obtained B+ FONLL reference is validated 
using published experimental cross sections measured in pp collisions at v^s = 7TeV Il20ll23l . 
The FONLL predictions are scaled by A(= 208), the atomic mass of the Pb nucleus, to account 
for the number of binary NN collisions Il42l . The FONLL uncertainties, which are larger than 
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the experimental uncertainties, represent the quadratic sum of several variations made to the 
calculation: of the factorization and renormalization scales, of the b quark mass, and of the 
uncertainty associated with PDFs (providing the largest contribution) Il37l - l3^ . The nuclear 
modification factor shown in Fig. 3 is computed as: 


nFONLL _ 

A ( WONTJ. ^ ' 

V dpx /pp 

where the numerator is defined in Eq. [^and the denominator is the corresponding theoretical 
calculation for B meson production in pp collisions at the same center-of-mass energy. The 
theoretical uncertainties represented by the open blue boxes in Fig. are computed by recal¬ 
culating with the upper and lower values of the FONLL predictions represented 

by dashed histograms in Fig. The nuclear modification factors of the three B mesons do not 
show evidence for modification of pPb data compared to the FONFF reference, in the consid¬ 
ered Pt range within the quoted uncertainties. No significant differences are observed between 
the three B meson species. In the lowest px interval measured, RpA^^^(pT) is 1.11 ± 0.08 (stat) ± 

0.17 (syst pPb) (syst FONFF), 0.99 ± 0.15 (stat) ± 0.18 (syst pPb) +° (syst FONFF), and 0.93 ± 
0.18 (stat) ± 0.20 (syst pPb) ('(g 24 (syst FONFF), for B+, B°, and Bg, respectively. 

The production cross section of B+ is also studied as a function of its rapidity in the center-of- 
mass frame (i/cm)- The ycM'differential cross section of B+ in the interval 10 < px < 60GeV/c 
is shown in Fig. (left). In Fig. (right), the rapidity dependence of the nuclear modification 
factor of B+ is shown. No strong evidence of rapidity dependence of Rp®^^^ is observed within 
the uncertainties. 

In summary, the first measurements of the B+, B'^, and Bg meson production cross sections in 
pPb collisions at = 5.02 TeV are presented. The mesons are measured in |yiab| < 2.4 and 
10 < Px < 60 GeV/c via the reconstruction of one of their exclusive hadronic decay channels. 
Within the transverse momentum and rapidity ranges studied, no significant modifications 
are observed, considering the statistical and systematical uncertainties, when compared to pp 
FONFF calculations scaled by the number of incoherent nucleon-nucleon collisions. These 
results provide a baseline for the study of in-medium b quark energy loss in PbPb collisions. 
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Figure 4: (left) The ycM'differential production cross section of B+ measured in pPb collisions 
at = 5.02 TeV. Vertical bars (the boxes) correspond to statistical (systematic) uncertain¬ 

ties. The listed global systematic uncertainty is not included in the data points. The result is 
compared to a FONLL calculation 113711391 represented by a continuous histogram. The dashed 
histograms represent the theoretical uncertainties for the FONLL reference, (right) The nuclear 
modification factor Rpa'^^(j/cm) of B+ as a function of ycM- The statistical and systematic 
uncertainties on the pPb data are shown as bars and yellow boxes around the data points, re¬ 
spectively. The systematic uncertainty from the FONLL reference is plotted separately as open 
blue boxes. The global systematic uncertainty is shown as a full grey box at unity, and is not 
included in the data points. 
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